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摘要 : 介绍 了 关于 暗物质 晕 的 边界 定义 问题 。 传 统 的 球形 高 密 算法 ， 将 一 定 密度 的 球形 高 密 
区 定义 为 一 个 上 暗 晕 ， 该 定义 方式 会 出 现 “ 伪 演化 ”的 问题 。 早 期 寻找 暗 晕 的 算法 还 包括 常用 的 连 
接 邻 近 粒 子 的 “朋友 链 ” 算 法 。 此 外 也 进行 过 各 种 探索 ， 如 : 根据 静态 区 间 的 最 远 位 置 而 定义 的 
音 坚 静态 半径 ， 以 及 近 些 年 根据 密度 轮 亡 在 外 围 剧烈 下 降 这 一 特点 而 提出 的 回 溅 半径 、 根 据 暗 坚 
的 动力 学 特征 所 提出 的 特征 耗竭 半径 和 内 部 耗竭 半径 ， 这 些 定义 能 够 反映 暗 坚 的 动态 特性 。 
还 介绍 了 计算 回 溅 半径 的 两 种 程序 和 观测 支持 ， 以 及 对 于 银河 系 的 内 部 耗竭 半径 的 探测 。 对 暗 坚 
的 并 合 树 以 及 增长 过 程 也 进行 了 介绍 。 最 后 对 各 种 边界 进行 总 结 和 展望 。 

x b dg: Wie, 边界 ， 高 密度 PRF: 耗竭 半径 
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1 引 言 


mu 


恨 据 目前 人 们 对 宇宙 大 尺度 结构 形成 的 理解 ， 宇 宙 是 由 大 约 27% 的 暗物质 、5%% 的 
子 物质 ， 以 及 其 余 68% 的 暗 能 量 组 成 上 。 冷 暗物质 分 级 拥 缩 冷 却 形成 暗物质 晕 (简称 上 暗 晕 )， 
而 瞳 党 是 组 成 宇宙 的 基石 ， 在 其 引力 势 的 中 心 孕 育 着 星系 。 

人 们 目前 对 于 暗 晕 的 基本 性 质 认 识 还 不 够 完善 ， 对 于 暗 学 的 大 小 这 一 重要 问题 的 认识 还 
处 在 起 步 阶段 。 最 初 ， 人 们 受到 球 对 称 南 缩 模型 ”的 启发 ， 提 出 了 球形 高 密 算法 (spherical 
overdensity, SO)””， 即 将 密度 达到 一 定 值 的 球形 区 域 定义 为 一 个 暗 量 ， 但 是 这 种 基于 密度 
的 定义 并 不 能 反映 暗 坚 的 真实 情况 ， 会 出 现 后 来 被 称 作 是 “ 伪 演 化 ”的 现象 ， 即 根据 定 
义 ， 瞳 晕 半 径 和 质量 增长 很 多 ， 而 实际 物理 过 程 并 没有 明显 变化 的 情况 。 此 外 ， 传 统 的 暗潮 
定义 方法 还 有 著名 的 “朋友 链 ”(Friends-of-Friends, FoF) 算法 中 下， 两 者 都 是 早期 定义 暗 
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后 来 ， 人 们 对 此 进行 了 不 少 探 索 ， 更 多 地 从 动态 的 视角 出 发 ， 而 不 是 简单 地 考虑 上 暗 坚 的 
EHE. Cuesta 等 人 (2008) ”提出 的 静态 半径 的 概念 ， 能 够 反映 暗 晕 的 动力 学 特性 ， 因 而 有 具 
定 的 优点 。 

前 些 年 ， 有 人 提出 暗物质 曼 的 回 溅 (splashback) 半径 (Rep) ONE mw 
Diemer 和 Kravtsov” 发 现 暗 晕 密 度 轮廓 的 斜率 在 外 围 处 剧烈 下 降 ， 该 半径 对 应 于 暗物质 
粒子 落 入 瞳 晕 后 第 一 次 经 过 轨道 远 心 点 的 位 置 。 瞳 晕 近 期 所 吸 积 的 物质 出 现 了 回 溅 ， 即 发 
溅 出 去 然后 又 回来 的 现象 。 该 半径 能 够 将 上 暗 曙 的 内 外 区 域 明 显 地 区 分 开 来 ， 因 而 被 认为 是 
非常 自然 和 物理 的 定义 ， 且 一 经 提出 就 备 受 关注 。 并 且 在 观测 上 也 能 找到 支持 ， 如 : More 
等 人 ”被 认为 是 最 早 测 得 回 溅 特征 的 ， 他 们 在 斯 隆 数字 巡天 (Sloan Digital Sky Survey, 
SDSS)" 中 星系 团 周围 的 测 光 星系 的 投影 数 密度 轮廓 中 发 现 了 一 个 尖锐 的 暗 晕 边缘 ， 而 
该 边缘 就 对 应 于 卫星 星系 落 入 后 的 第 一 轨道 远 心 点 。 此 外 ，Chang 等 人 (2018)” 运 用 暗 能 
量 巡 天 (Dark Energy Survey, DES) T 的 数据 ， 测 得 了 星系 密度 轮廓 的 类 回 溅 的 陡峭 的 特 
iE, Bianconi 等 人 (2021) ^ 利用 近邻 星系 团子 结构 巡天 (Local Cluster Substructure Survey, 
LoCuSs)™ 的 数据 ， 直 接 探测 了 星系 团 的 回 溅 特征 ， 并 研究 了 星系 团 的 动力 学 状态 对 该 特 
征 出 现 的 影响 。 

无 独 有 偶 ， 最 近 Fong 和 Han” 引入 了 暗 晕 的 耗竭 半径 (depletion radius) 的 概念 ， 分 为 

内 部 耗竭 半径 (inner depletion radius, Ria) 和 特征 耗竭 半径 (characteristic depletion radius, 
.a)， 其 定义 分 别 是 在 质量 流入 率 最 大 的 位 置 和 偏 祖 轮 廓 最 小 的 位 置 。 其 中 的 内 部 耗竭 半径 
(Ria) 将 正在 增长 中 的 暗 蛙 ， 与 周围 暗 军 因 吸 积 而 被 “消耗 ”的 环境 分 隔 开 来 ， 并 且 该 区 域 
内 的 密度 具有 大 致 相同 的 值 。 随 后 ，Li 和 Han (2021)” 利 用 近邻 矮星 系 的 径 向 运动 ， 率 先 
测 出 银河 系 暗 尝 的 内 部 耗竭 半径 和 回转 半径 (turnaround radius, Rta)» Hf HABA T 4R 
河 系 的 独特 演化 历史 和 动力 学 状态 。 
第 2 章 将 介绍 几 种 关于 暗 军 边界 定义 的 方法 ， 尤 其 是 传统 的 球形 高 密 算 法 及 其 问题 ， 并 
重点 介绍 回 溅 半径 的 发 现 过 程 、 计 算 回 溅 半径 的 程序 算法 以 及 耗竭 半径 的 引入 过 程 ; 第 3 
章 将 介绍 一 些 观测 上 用 到 的 和 观测 到 的 瞳 军 边界 ; 第 4 章 将 介绍 瞳 晕 的 并 合 树 (merger tree) 
以 及 对 于 暗 晕 增 长 过 程 的 描述 ; 第 5 章 将 对 各 种 边界 进行 总 结 和 探讨 。 
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O 32 暗物质 晕 边 界 的 几 种 定义 


本 章 主要 介绍 暗 党 边界 的 几 种 主要 的 定义 ， 分 别 是 传统 的 球形 高 密 算法 (SO) 和 FoF 算 
法 ， 以 及 后 来 提出 的 静态 半径 、 回 溅 半径 和 耗竭 半径 。 基 于 密度 的 算法 存在 “ 伪 演 化 ”的 
问题 ， 难 以 反映 上 暗 晕 的 真实 情况 。 而 著名 的 FoF 算法 ， 通 过 连接 数值 模拟 中 一 定 距 离 以 内 
的 粒子 ， 以 此 来 定义 瞳 晕 ， 并 且 ， 很 多 追踪 上 暗 晕 和 构建 并 合 树 的 程序 往往 先 运 行 FoF 算法 。 
此 外 ， 人 们 也 试图 寻找 新 的 办 法 去 定义 瞳 坚 的 边界 。 根 据 径 向 速度 提出 的 静态 半径 和 相应 的 
静态 质量 的 概念 ， 也 具有 一 定 的 优点 。 本 章 重 点 介绍 回 溅 半径 的 概念 ， 人 们 发 现在 暗 晕 的 外 
围 密度 轮廓 的 (对 数 ) 斜率 是 剧烈 下 降 的， 经 过 深入 研究 而 定义 出 了 这 一 著名 的 概念 。 我 们 
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也 将 展开 介绍 两 种 计算 回 溅 半径 的 程序 。 最 近 又 有 人 根据 暗 晕 的 偏 祖 轮 廊 和 物质 在 径 向 流 
动 的 动态 过 程 ， 定 义 了 特征 耗竭 半径 和 内 部 耗竭 半径 这 一 组 物理 量 ， 并 且 认 为 这 更 加 自然 地 
刻画 了 暗 晕 边 界 ， 比 回 溅 半径 意义 更 为 清晰 。 此 外 ， 需 要 强调 的 是 ， 伪 演化 问题 是 密度 定义 
所 独 有 的 ， 其 他 反映 瞳 晕 动态 特征 的 定义 则 不 存在 此 问题 。 
2.1 ”基于 密度 的 定义 
在 球 对 称 南 缩 模型 "的 启发 下 ， 人 们 发 明了 球形 高 密 算法 (SO)””， 即 将 密度 达到 一 定 
值 的 球 对 称 区 域 当做 一 个 上 暗 晕 。 
具体 过 程 为 ，SO 算法 先 将 数值 模拟 中 的 所 有 粒子 根据 它们 的 局 部 密度 排序 ， 然 后 从 密 
度 最 高 的 地 方 出 发 ， 划 出 球形 的 区 域 ， 该 区 域 半径 不 断 增长 ， 直 到 其 密度 下 降 到 一 定 的 阔 值 
以 下 。 将 这 团 暗 党 物质 的 中 心 根据 所 选区 域 的 质量 中 心 重新 运 代 ， 然 后 重新 划 出 球形 区 域 ， 
直到 上 暗 曼 的 中 心 和 成 员 收敛 为 止 。 随 后 ， 所 有 已 经 被 分 配 到 暗 晕 的 粒子 被 剔除 出 去 不 再 计 
算 。 用 同样 的 方法 反复 欠 代 ， 直 到 找 出 所 有 上 暗 棠 。 最 终 ， 内 部 小 的 暗 晕 或 者 与 更 大 的 暗 曼 重 
县 的 暗 晕 将 合并 成 为 更 大 的 团 块 。 
通过 以 上 方法 找到 的 上 暗 晕 ， 其 质量 、 密 度 和 半径 符合 以 下 公式 ; 


4n 
Ma = Ane RÀ , (1) 


FOF, prop 叫做 参考 密度 ，A 则 是 倍数 。 参 考 密 度 可 取 宇 宙 物 质 平均 密度 或 者 临界 密度 ， 
这 两 种 密度 都 会 随 红 移 演化 。 根 据 选择 有 200pmean 和 180puit 等 ， 相 应 的 半径 为 Room 和 
Rigoe; 质量 为 M»90m 和 Mhi0ce 和 常用 的 是 用 Bryan 和 Norman"? 1998 年 提出 的 密度 (简称 


BN98): 
pvir = (187? = 82q = 39q”) Perit ; (2) 
其 中 ， 
(2 
q = = (3) 


Nm oa73 + Nao 
与 之 相对 应 的 就 是 Rs 和 Mi;.， 即 人 们 通常 所 说 的 位 力 半径 和 位 力 质 量 。 式 中 a 是 尺度 因 
d^ pau EMERE, Ruo 和 ao 分 别 是 红 移 为 0 时 的 宇宙 物质 密度 参数 和 宇宙 学 常数 
参数 。 

然而 ，SO 算法 试图 通过 密度 区 域 来 定义 暗 曙 边界 存在 着 一 定 的 问题 ， 不 能 够 反映 暗 尝 
真实 过 程 ， 会 出 现 所 谓 的 “ 伪 演 化 ”的 现象 ， 即 根据 定义 ， 暗 晕 的 质量 和 半径 增长 很 多 ， 而 
实际 并 未 伴随 内 部 的 物理 过 程 也 变化 的 现象 。 以 下 将 详细 介绍 伪 演 化 现象 。 

Diemand 等 人 ”于 2007 年 给 出 在 Via Lactea 数值 模拟 ”中 银河 系 质量 的 暗 晕 质 量 增 
KA, WAG 从 图 四 中 可 以 看 出 ， 物 理 质量 的 吸 积 过 程 从 大 约 红 移 z = 1.7 (尺度 因子 
a = 0.37) 以 后 就 很 少 发 生 ， 并 且 在 400 kpc 以 内 的 尺度 范围 内 ， 红 移 z — 0 时 80% 以 上 的 
物质 在 红 移 z = 1 的 时 候 就 已 经 吸 积 完成 。 而 在 相应 的 时 间 范 围 里 ， 其 位 力 质量 My 以 及 
JMz00 却 增 长 了 很 多 。 

Cuesta 等 人 "给 出 星系 尺度 的 瞳 晕 密度 轮廓 演化 OLA). MAB 中 可 以 清晰 地 看 出 ， 
自从 红 移 z = 1 以 来 ， 暗 晕 的 密度 轮廓 几乎 没什么 变化 ， 但 是 红 移 z = 0 时 的 位 力 半径 却 增 
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M(«r)/10?M, 


oC ENEE EMM J 
0.2 0.4 0.6 0.8 1 
a=1/(1+z) 
VE: 上 方 红色 虚线 为 Maoo， 下 方 红色 虚线 为 Mi;:， 蓝 色 点 线 为 最 大 旋转 速度 半径 内 的 质量 ， 一 条 条 黑色 的 实 线 


为 一 定 物理 半径 的 球形 区 域内 的 质量 演化 。 


1 ， 数值 模拟 中 一 个 银河 系 质量 大 小 的 暗物质 晕 的 质量 吸 积 历史 “ 


加 到 了 > = 1 时 的 两 倍 。 


108 


10°} 


104 


1012 


p(r)/(h’- Mo Mpc”) 


r/(h?-Mpc) | 


VE: 竖 直 线 是 各 个 红 移 时 的 位 力 半径 。 


图 2 ”数值 模拟 中 一 个 星系 大 小 的 暗物质 晕 的 密度 轮廓 的 演化 ” 


对 此 ，Diemer 等 人 ”研究 发 现 ， 红 移 z = 1 以 后 ， 相 当 多 的 暗 晕 存 在 伪 演 化 现象 ， 并 
且 大 部 分 小 质量 (< 10 Mo/h) 的 暗 晕 都 有 伪 演 化 的 发 生 。 
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2.2 “朋友 链 ”(FoF) 算法 


FoF 算法 ”将 距离 小 于 一 定 值 on-13 的 粒子 连接 起 来 ， 其 中 n 是 平均 粒子 密度 ， 通 


ChinaXiv 合 作 期 


539 


HAL b = 0.2， 即 粒子 平均 距离 的 20%。 由 此 便 得 到 一 条 条 粒子 的 链条 ， 将 足够 质量 的 作为 


Hie. FoF 算法 的 一 个 特征 是 ， 计 算得 到 的 瞳 晕 形状 不 规则 ， 但 是 偶尔 会 


接 的 粒子 将 两 个 暗 坚 连 接 起 来 的 情况 。 当 然 ， 该 方法 只 是 定义 了 暗 坚 本 身 
“Re”. 


FoF SE 72 18 Bx ze A EHA W SIS A | WSS 


SUBFIND)””、 等 级 型 束缚 追踪 算法 (Hierarchical Bound-Tracing, HBT) 


出 现 一 些 中 间 连 
， 不 存在 严格 的 


(subhalo finder, 


EAEX K 空间 拓 


扑 适 应 改进 的 稳健 高 密度 算法 (Robust Overdensity Calculation using K-Space Topologically 


Adaptive Refinement, ROCKSTARJ) 呈 、 子 晕 链 接 算法 (SUBLINK) 等 。 
2.3 ”静态 半径 


Cuesta 等 人 ”给 出 数值 模拟 中 三 个 质量 区 间 的 暗 晕 平均 径 向 速度 ， 各 


区 间 由 数 百 个 不 


同 的 暗 坚 取 平均 得 到 ( 见 图 团 。 上 下 虚线 之 间 是 他 们 所 划分 的 静态 区 间 (5% 的 位 力 速度 )。 


iv 


Fe AR Hr tt FE Hid 22 HUE TRE RII A SR ( 径 向 速度 为 负 )， 强 内 落 以 外 的 区 域 ， 即 使 平均 径 


向 速度 接近 零 ， 也 不 属于 静态 区 间 。 小 质量 和 星系 质量 的 暗 晕 则 在 哈 勃 流 之 前 显示 出 小 的 外 


流 。 他 们 将 静态 区 间 的 最 远 位 置 定义 为 静态 半径 ， 其 中 所 包含 的 质量 即 为 静态 质量 ， 并 认为 


这 比 基 于 SO 算法 的 位 力 质量 和 半径 更 好 。 


1 T T —— ty T T a a AE 


—— 10" h-M,<M,,,<2x10" ht Mg 
0.8[- — — 10? i'- Me «M, «2x10? I- M 
----104 bi Me«M,, 2x104 Ir Me 


0.6 


vir 


0.4 


<V>/V， 


ik: 短 虚 线 为 星系 团 质 量 大 小 的 暗 晕 ， 长 虚线 为 星系 质量 大 小 的 暗 晕 ， 实 线 为 小 质量 暗 曙 。 


图 3 不同 质量 区 间 的 暗物质 景 的 平均 径 向 速度 中 


2.4 BREE 


前 几 年 提出 的 回 溅 半径 (Rso) 的 概念 ”” ， 被 认为 是 较为 自然 的 暗 晕 边 界 。 
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Tal yA 
ZH 


Diemer 和 Kravtsov 出 红 移 z = 0 时 较 大 质量 的 暗 虹 (v > 3.5, Mor > HMM) 
S PE FER FHE I T B ( 见 图 团 。 可 以 看 出 该 轮 廊 在 外 围 处 剧烈 下 降 ， 与 传统 的 NEW 和 
Einasto” 轮廓 对 比 相 去 其 远 。 


dlgp/dlgr 


y= 


0.1 05 1 5 


VE: 实 线 是 较 大 质量 暗 晕 密 度 轮 廓 的 斜率 ， 点 划 线 和 虚线 分 别 是 NFW 和 Einasto 轮廓 。 
4 BERE SEU 


对 于 斜率 跳跃 下 降 的 原因 ，Adhikari 等 人 ”总 结 道 ， 局 部 陡峭 是 由 壳 层 交叉 形成 的 焦 
散 线 所 导致 ， 与 暗 晕 近期 所 吸 积 物质 的 回 溅 现象 有 关 ( 即 飞溅 出 去 然后 又 回来 ， 如 同 防 溅 挡 
板 一 样 )。 各 个 粒子 的 轨道 在 相同 的 地 方 ， 尤 其 是 轨道 远 心 点 附近 产生 堆积 ， 使 得 焦 散 线 加 
强 。 另 外 ， 这 一 显著 的 焦 散 线 在 球 对 称 ” 和 非 对 称 ” 的 情况 下 都 存在 。 据 此 ， 人 们 便 定义 
了 暗 坚 的 回 溅 半径 。 

More 等 人 中 给 出 红 移 > = 0 时 ， 质量 Mi = 11x104Mgo/h,. REMAX = 0.8 
的 上 暗 雷 的 切片 投影 密度 ( 见 图 团 。 从 图 加 中 可 以 清晰 地 看 出 ， 密 度 场 在 Rp 的 地 方 确实 是 
剧烈 下 降 的 。 这 样 的 回 溅 半径 将 暗 晕 的 多 流 区 域 与 内 落 区 域 有 效 地 区 分 开 来 。 因 此 ，R,, 被 
认为 是 暗 晕 边界 的 一 个 非常 自然 的 定义 。 同 时 ，Rse 以 内 的 质量 就 相应 地 定义 为 回 溅 质量 ， 
Msp = M(« Rs). 

More 5& A" ggg zs EOE AU ASE ABE A TE TE LE. fi HB HT ASE 
回 溅 半径 的 拟 合 公式 : 

Rep 
R200m 


ABD 是 质量 吸 积 率 ， 表 示 为 : 


= 0.54[1 + 0.53(2,(2)](1 +1.36e-7/3.04) , (4) 


I = Alg Mi/Alga , (5) 


Rooom 是 式 (M) 中 取 200 倍 宇宙 物质 平均 密度 定义 的 半径 ，a ERENT, Mair 是 位 力 质 
i, Q(z) 是 宇宙 物质 密度 参数 。 

此 外 ，Diemer 等 人 衬 守 开发 了 一 个 计算 回 溅 半径 的 程序 ， 子 党 和 粒子 轨迹 分 析 (Sub- 
halo and PARticle Trajectory Analysis, SPARTA) 算法 。 该 算法 分 析 了 宇宙 学 模拟 中 宿主 


[ Samet 


202306.00403v1 


chinaXiv 


4 期 


ik: 图 中 白色 的 实 线 为 Ra. 


图 5 


章 坚 中 的 全 部 粒子 轨道 ， 并 且 能 准确 提取 出 粒子 落 入 暗 晕 后 的 第 一 远 心 点 ， 即 回 溅 的 位 置 。 
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lgp/p。 


X/( Mpc) 
点 划 线 为 R200m> 虚线 为 nm 点 线 为 Rintatie 
BEA 0.15Ro0m 的 瞳 物质 量 切 片 的 投影 密度 于 


他 们 将 粒子 的 远 心 点 半径 的 平滑 平均 值 定义 为 暗 尝 的 回 溅 半径 。 


图 日 给 出 数值 模拟 中 


图 日 展示 了 传统 的 位 力 半 径 和 用 SPARTA 算法 给 出 的 回 溅 半径 。 
子 数 Nooom > 1000， 相 应 的 质量 在 1.4 x 10 Mo/h 及 以 上 。 其 中 


M200m =1.2~x 104 Mo /he 


的 一 个 宽度 和 深度 均 为 30 Mpc/h， 高 度 为 15 Mpc/h 的 切片 。 
Efe. Fr thon mm wm du 
， 中 间 大 的 暗 晕 质量 


注 ， 白 色 圆 为 Re WAGE 


图 6 


为 Ras 


SPARTA 算法 给 出 的 回 溅 半径 和 位 力 半径 对 比 


另外 ，Mansfield AT FR TRER ZI (SHELL Finding in Spheroidal Halos, 
SHELLFISH) 算法 。 该 算法 可 以 根据 数值 模拟 中 密度 场 的 单个 快照 ， 来 识别 出 暗 晕 的 “ 回 
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溅 壳 层 ”， 无 需 粒 子 的 动力 学 信息 。 
图 四 展示 了 SHELLFISH 算法 运行 的 4 个 步骤 。 


10° 
102 
M | 
Q 
= 
a 10 
107 
10°93 10 20 
R/(b^-Mpc) 
b) 


d) 


ik: SHELLFISH 的 步 又 。 


图 7 “对 一 个 星系 团 大 小 的 暗 晕 使 用 SHELLFISH 算法 加 


(1) 从 中 心 出 发 ， 随 机 作出 几 万 条 射线 (图 中 红色 的 )， 来 追踪 该 暗 晕 的 密度 场 。 
(2) 测量 此 射线 方向 上 的 密度 轮廓 ， 得 到 图 中 黑色 的 线 ， 经 过 Savitzky-Golay 过 滤 
器 四 平滑 后 ， 得 到 图 中 红色 的 线 ， 图 中 的 箭头 指向 平滑 轮廓 最 陡 的 地 方 。 
(3) 标 出 所 有 斜率 最 陡 的 点 ， 并 将 其 分 类 ;白色 曲线 是 在 点 的 选择 过 程 中 产生 的 过 滤 样 
条 (filtering spline); 距离 该 曲线 近 的 点 标 为 白色 ， 远 的 标 为 红色 。 
Ta (4) 根据 100 个 随机 方向 的 平面 中 总 体 回 溅 点 的 分 布 ， 拟 合 出 Penna-Dines 曲面 中 。 
Bi 由 此 SHELLFISH (875 3) T mri ( lex, dbi fc JIU RT DA 5 8 eL S 


Has ARR " 
Py = (Œ) (6) 
KB. Vp 为 壳 层 所 包含 的 体积 
2.5 ”耗竭 半径 (depletion radius) 
最 近 ，Fong 和 Han” peh T rs Pg BB REMAKE (inner depletion radius, Ria) 和 特征 
耗竭 半径 (characteristic depletion radius, Req) 这 一 组 全 新 的 概念 ， 并 且 也 认为 是 对 于 暗 晕 
边界 的 自然 刻画 


VE 
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图 四 的 纵 坐 标 b 为 偏 祖 因子 ， 表 示 为 ; 
&e(r) | (p) 

= 7 
Emm) Enr) ， (0) 
HFP, fam 和 £s 是 瞳 晕 -物质 、 物 质 -物质 相关 函数 ，6(r) 是 暗 晕 周围 的 物质 密度 差 。 图 中 
清晰 可 见 的 是 中 间 一 个 沟 岁 的 存在 ， 据 此 ， 他 们 便 定 义 了 暗 晕 的 特征 耗 剖 半 径 (Rea), Ba 
HAT b 最 小 的 地 方 。 


b(r) = 


a. 


10° 10! 10° 101 10° 10! 
r/(h?^-Mpc) 


注 : 不 同 的 子 图 代表 不 同 的 参量 ， 不 同 颜色 表示 不 同 的 质量 区 间 。 
图 8 ， 瞳 晕 的 偏 祖 轮 廊 与 半径 的 关系 呈 


图 回 展示 了 平均 偏 祖 因子 、 径 向 速度 和 质量 流速 (mass flow rate, MFR) 轮廓 ， 给 出 上 暗 
尝 的 质量 区 间 和 相应 区 间 内 的 瞳 晕 数目 。 全 部 的 径 向 速度 可 以 分 解 为 办 = wp ton, SSA 
边 的 两 项 分 别 是 本 动 速度 和 哈 勃 速度 ， 速 度 为 负 则 表示 朝向 中 心 。MEFR 表示 单位 时 间 流 
经 横 截 面 的 物质 ，MFR = p(ryu(r)4nr?, FA Mar: vvir TG. 进行 归 一 化 ， 其 中 位 力 速 度 为 
wir = V GMyara G 是 引力 常数 。 由 图 日 中 明显 的 沟 蜜 可见， 大 质量 上 暗 蜡 活跃 地 吸 积 
周围 物质 ， 而 小 质量 的 则 不 然 。 
由 于 MPR 决定 着 密度 轮廓 的 演化 ， 因 此 他 们 把 图 四 所 示 的 最 大 内 流 位 置 叫做 最 大 内 
流 半 径 。 在 最 大 内 流 半径 内 ， 物 质 像 是 被 倒 入 上 暗 举 ; 而 在 该 半径 外 ， 物 质 像 是 被 渐渐 抽空 一 
般 。 最 大 内 流 位 置 象征 着 从 暗 量 被 构建 的 过 程 到 物质 被 吸 积 所 消耗 的 过 程 的 转变 。 可 见 ， 最 
ACE i. 暗 晕 周围 活跃 的 消耗 区 的 内 边界 。 于 是 ， 偏 祖 因 子 最 小 处 的 特征 耗 曾 半 
， 应 当 在 最 大 内 流 半径 之 外 。 正 如 图 日 所 示 ，Rca 比 最 大 内 流 半 径 大 了 10% ~ 2096. d& 
is EIE, dd dE — dg X P CAE fe nu (lo “TA BRE ABE” (Ria) Ria 预示 着 连绵 
不 绝 的 耗竭 过 程 的 开始 ， 而 最 小 偏 祖 位 置 的 Ru 则 像 是 特征 或 者 “最 深 的 耗竭 ”半径 一 样 ， 
两 者 共同 刻画 了 耗竭 沟 蜜 的 尺度 大 小 。 
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& 
Q 
M 
= —— [11.50,12.05]:1 356 103 
- —— [12.05,12.60]:445 201 
m" -0.2| . [12.60,13.15].140 759 
(13.15,13.70]:41 114 
-0.3|  [13.70,14.25]:10 077 \ J/ 
(14.25,14.80]:1 640 \ 
[14.80,15.35]:104 j 
-0.4 
- 0.00 
$ 
> 
= -0.04 
= 
s 
E 
—0.08 
E 
= 
-0.12 
10? 109 
r/(h?^-Mpc) 


10! 


: 上 图 为 图 加 中 的 偏 祖 轮廓 ， 中 图 和 下 图 分 别 是 径 向 速度 和 MPR. Elm np 
十 字符 号 代表 特征 耗竭 半径 和 
不 同位 力 质量 区 间 的 暗 晕 平 均 偏 祖 因 子 、 径 向 速度 和 M FR 轮廓 图 


Hu 


HEME, 


色 代 表 质 量 区 间 ， 五 角 星 和 


— t 


0 
一 一 tt 


消耗 


: Sa 
<+<+<+ «——— 4—— 4— - 


> 


r 


T 


的 物质 流动 的 幅度 。 


VE: 从 黑 实 线 到 红 虚 线 代 表 了 暗 晕 的 演化 ， 垂 直 点 线 表 示 最 大 内 流 的 位 置 即 内 部 耗 章 边 界 ， 水 平 箭头 代表 最 大 内 
流 位 置 周围 


图 10 。 暗 晕 在 最 大 内 流 位 置 附近 吸 积 的 示意 图 加 
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Fong 和 Han 给 出 不 同 质量 区 间 的 暗 晕 的 径 向 速度 分 布 ( 见 图 J)， 位 力 质量 Mv 以 


Mg -h-! 为 单位 ， 在 他 们 此 项 工作 中 ， 


度 轮廓 的 相 空间 分 布 最 窄 的 位 置 。 


音 尝 所 在 红 移 均 设 为 z = 0。 图 中 的 耗竭 半径 位 于 密 


4 — 
[11.50,12.05] e 
LI al, LI 


VD/ Voi 


05 15 3235 35 05 


a) 
YE: a), b), c) 表示 以 暗 晕 的 位 力 质量 取 对 数 (lg 


ti I 了 | 
15 25 35 05 15 25 3.5 
lg[r/(kpc-h*)] 
b) 


Muvi) 划 分 的 暗 晕 区 间 ， 图 中 只 显示 每 个 区 间 随 机 选择 的 100 


个 瞳 晕 。 颜 色 图 表示 粒子 的 分 布 ， 由 外 到 内 的 白色 


线 包含 99%, 80%, 60% 的 粒子 的 等 密度 


ae 
E 
3E 

B 
l2 
ou 
St 

B 


线 分 别 是 径 向 速度 和 本 动 速 度 ， 垂 直 的 红 实 线 、 


品 红 色 点 划 线 、 蓝 色 虚 线 分 别 是 耗竭 半径 、 回 溅 半径 、 位 力 半 


径 的 位 置 ， 蓝 绿色 的 又 号 代表 回转 半径 的 位 置 ， 黑 


色 点 标记 最 大 内 落 的 位 置 。 


图 11 不同 质量 区 间 的 归 一 化 径 向 速度 分 布 吧 


BLEW REM A> Ria 


FY DA EL Sls AOA ee TL AE TR B] doe Vb A OP B e E fe. 
便 在 距离 最 窄 速 度 分 布 边界 不 远 ， 但 更 内 部 的 位 置 。 


为 了 更 清晰 地 说 明 ， 他 们 在 图 四 中 给 出 了 一 个 Me = 1.35 x 10 Mo - 7 FESTE 


向 速度 相 空 间 分 布 ， 


可 见 最 大 内 落 的 位 置 清 晰 地 标示 出 回 溅 边界 的 外 围 。 


所 以 ， 内 部 耗竭 边 


界 可 以 理解 为 包含 了 完整 的 回 溅 族群 的 边界 。 


3 ”上 暗物质 尝 几 种 边界 的 观测 研究 和 应 用 


本 章 将 介绍 观测 上 常用 的 边界 Rsoo， 


以 及 关于 回 溅 半径 (Rep) 和 内 部 耗竭 半径 (Ria) 的 


观测 支持 。 对 于 暗 晕 边界 回 溅 半径 的 密度 轮廓 剧烈 下 降 的 现象 ， 观 测 上 ， 已 经 发 现 星系 团 中 


星系 的 数 密度 轮廓 和 物质 质量 轮廓 也 存在 类 似 现 象 ， 而 在 耗竭 半径 这 一 


河 系 的 内 部 耗竭 半径 也 有 人 将 其 具体 测 出 
3.1 ”基于 密度 定义 的 Rsoo 的 观测 应 用 


观测 上 ， 人 们 根据 苏 尼 亚 耶 夫 - 泽 尔 多维 奇 (Sunyaev-Zel'dovich, SZ) 效应 “估计 质量 
时 ， 经 常用 Mso 和 Rsoo 来 衡量 星系 团 的 质量 和 大 小 。 


概念 提出 后 不 久 ， 银 


如 Ricci 等 人 (2020)"" H XXL iW 
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ik: 图 中 灰色 的 点 表示 随机 选取 瞳 晕 中 1% 的 粒子 ， 黑 色 和 灰色 曲线 是 全 部 的 和 本 动 径 向 速度 的 平均 值 ， 两 条 垂 
直线 标记 的 是 位 力 半 径 和 最 大 内 落 速度 的 位 置 。 


12 ”单个 瞳 晕 的 径 向 速度 分 布 图 加 


XU 的 数据 研究 XLSSC 102 星系 团 呈 时 ， 就 使 用 了 式 四) 的 关系 式 ， 取 500pwrit， 即 该 星系 
团 所 在 红 移 处 500 倍 的 字 宙 临界 密度 ， 如 图 T3 所 示 。 


ee 


ve: 图 中 黑色 虚线 圆圈 是 Rsoo， 背 景 是 射线 面 亮度 图 ， 绿 色 的 是 SZ 信 品 比 (S/N) 轮廓 线 ， 品 红色 的 线 是 
oo (S/B) 轮廓 线 ， 品 红色 的 又 号 (x) 是 最 亮 团 星 系 (BCG) 的 位 置 ， 粉 红色 的 箭头 代表 可 
4 并 合 轴 。 


T 


13 ”星系 团 XLSSC 102 的 多 波段 图 像 中 
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3.2 [IEEE (Ro) 的 观测 支持 

2016 年 ，More 4& A" x Jc gu fe i e 
的 回 溅 半径 。 他 们 利用 SDSS DR8” 的 星 
系 团 周 围 星 系 的 面 密度 轮廓 (X, (R)), È 
现 类 似 于 Diemer 和 Kravtsov £8 $8] hy sb 
率 剧烈 下 降 的 现象 ， 如 图 四 所 示 。 对 于 图 
中 观测 与 模拟 中 回 溅 半径 的 差异 ， 一 方面 
可 能 因为 参数 测量 偏差 ， 另 一 方面 也 可 能 
是 由 于 暗物质 粒子 的 自 相互 作用 包含 着 未 
知 的 新 物理 ， 从 而 导致 回 溅 半径 的 减 小 。 
此 外 ，Chang 等 人 (2018) "iz DES 
的 数据 ， 研 究 了 星系 团 周围 的 星系 数 密度 BE i) 
和 弱 引 力 透 镜 质量 轮廓 ， 也 得 到 了 星系 密 biis ei 
FE t: ECKL LIRA EE, um ERE 中 14 DE US OSPITI SESS 
所 示 。 而 且 ， 重 要 的 是 ， 他 们 的 弱 引力 透 ar 
镜 测 量 ， 第 一 次 测 得 星系 团 物质 轮廓 的 这 一 现象 。 对 于 图 吧 中 星系 数 密度 的 轮廓 和 弱 引 
力 透 镜 得 到 的 物质 质量 轮廓 ， 斜 率 最 陡 的 位 置 相互 符合 ， 两 者 的 幅度 和 形状 也 一 致 ， 均 比 
NEW 轮廓 要 陡峭 ， 符 合 对 于 回 溅 半径 特征 的 描述 。 


dlgz /dlgR 


& -2 

"UO 
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x 

E - 

3 -3 
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ik: 图 中 灰色 的 带 是 星系 的 数 密度 轮廓 ， 红 色 的 带 是 弱 引 力 透 镜 得 到 的 质量 轮廓 ， 垂 直线 是 各 自 的 平均 Rep» 7K 
平 误差 棒 是 各 自 的 不 确定 度 。 


图 15 ”星系 团 的 星系 数 密度 和 物质 质量 轮廓 的 斜率 图 


3.3 ”内 部 耗竭 半径 (Ria) 的 观测 研究 

2021 年 ， 在 内 部 耗竭 半径 (Ria) 提出 后 不 久 ，Li 和 Han 便利 用 银河 系 附近 3 Mpc 内 
itg ” 和 附近 的 矮星 系 ” 的 目录 ， 并 根据 附近 可 用 的 矮星 系 ” 的 自行 运动 ， 将 相对 于 
日 心 的 视 向 速度 转换 为 相对 于 银 心 的 径 向 速度 (不 可 用 的 则 忽略 自行 )， 首 次 测 出 银河 系 暗 
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SERA BEES e (Ria)。 男 外 ， 他 们 顺便 也 测定 了 其 回转 半径 (Ria), BISERI A E E 
达到 的 径 向 速度 为 零 的 位 置 。 

图 上 中 ， 与 M31 的 距离 在 600 kpe 以 内 的 星系 用 空心 圆 标 出 ， 在 分 析 时 忽略 。 他 
们 用 其 余 的 星系 (实心 圆 ) 计算 了 平均 速度 轮廓 (绿色 的 实 线 ) 和 相应 的 lc 的 不 确定 度 
(绿色 的 带 )。 最 终 ， 他 们 得 到 的 银河 系 内 部 耗竭 半径 Ru = (559+ 107) kpc, RABE 
Ria = (839 + 121) kpc. 


位 力 化 增长 ”消耗 后 退 
© 
© 
— 100+ o Oo 
oO 
is p. 8 © o 9 ony LI 
g © CERE & Q O @ 
EJ Um m c up ode 1. 
© © 让， 
f © 0。 eo% Oe? 
-100 自 © [3] O [e] HM31 o 
$ o8» 
-200r @ 6 Ó 
© 
i. L |: V-Hr 
aon po a a acad : 1 PE © hana am a 
30 100 300 1 000 3 000 
r/kpc 


ik: 图 中 的 五 角 星 代 表 测 得 的 内 部 耗竭 半径 (Ria) 和 回转 半径 (Ri). 
图 16 ”银河 系 附 近 3 Mpc 内 的 星系 的 径 向 速度 分 布 四 


4 暗物质 晕 的 演化 过 程 


本 章 介绍 反映 暗 晕 并 合 过 程 的 并 合 树 ， 这 是 理解 暗 晕 结构 绕 不 开 的 话题 。 还 举例 介绍 位 
力 半径 、FoF 暗 晕 、 静 态 半 径 和 回 溅 半径 所 定义 的 暗 曙 增长。 位 力 半 径 的 情况 下 可 以 明显 
看 出 暗 尝 的 增长 过 程 ， 分 为 早期 快速 增长 和 晚期 慢 速 增长 两 个 阶段 。 其 他 几 种 定义 也 各 有 特 
点 ， 一 同 列 出 以 便 参考 。 
4.1 暗物质 晤 的 并 合 树 

数值 模拟 表明 ， 当 两 个 暗 曼 并 合 时 ， 小 质量 的 暗 曼 并 不 会 完全 消失 ， 而 是 会 继续 以 一 
种 自 束缚 的 ， 叫 做 子 曙 (subhalo) 的 形态 存在 ， 子 曙 受 到 大 质量 暗 曙 引力 势 的 作用 而 绕 其 旋 
te, AIO 展示 了 暗 晕 的 一 连 串 的 并 合 过 程 ， 紫 色 的 球体 是 主 祖先 (main progenitor) ÆW 
并 合 历史 ， 称 之 为 零 级 祖先 (zeroth-order progenitor)。 零 级 祖先 吸 积 一 级 祖先 (first-order 
progenitor) (图 中 检 色 的 球体 )， 后 者 在 红 移 > = zo 时 变 成 一 级 子 晕 (first-order subhalo). 


mE 
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+a 


ik: FER PRA CE UGE MF Ae, REO. 
记 这 是 一 次 平滑 吸 积 。 


色 分 别 代表 一 级 、 二 级 祖先 ( 子 景 )， 小 的 灰色 方块 标 


Er: 


17 —ABES (host halo) 的 并 合 树 示意 图 加 


同 理 ， 一 级 祖先 吸 积 二 级 祖先 (second-order progenitor) (KH EA 中 蓝 色 的 球体 )， 使 其 在 红 移 
z = zo PI AE, 2A T (second-order subhalo) 或 者 叫 子 子 军 (sub-subhalo)。 以 此 类 推 ， 追 
踪 主 祖先 晕 的 分 支 叫做 主 分 支 (main branch)”. ifj 2 节 提 到 的 几 种 瞳 晕 追踪 和 构建 并 合 
树 的 算法 ， 首 先 通过 FoF 算法 找到 暗物质 团 块 ， 然 后 根据 各 自 的 标准 ， 得 到 暗 蛇 和 图 吧 所 
ZR HERE EAE APT. 
4.2 ”暗物质 晕 的 增长 过 程 
4.2.1 位 力 半 径 的 描述 

图 四 是 Zhao 等 人 ” 给 出 的 暗 晕 半 径 的 增长 过 程 。 从 图 中 可 见 ， 三 种 不 同 的 密度 定义 
的 半径 差别 不 大 。 暗 党 常用 的 NFW 轮廓 ”中 ， 内 部 特征 半径 rs 定义 在 对 数 密度 斜率 为 -2 
的 位 置 ， 聚 集 度 (concentration) c = ryir/Ts rv 为 位 力 半径 。 

图 四 是 Zhao 等 人 ”给 出 的 主 分 支 上 的 暗 晕 质量 增长 历史 。 对 于 瞳 晕 的 平均 增长 历史 ， 
他 们 给 出 如 下 公式 : 


Myir(2 aps 

ee I apar te ^ (8) 

HP, £= pvirtp/Ppviz(z)， 以 此 来 表示 “时 间 ”， 下 标 tp 表示 转折 点 (turning point), pvir 
表示 位 力 密度 ， 对 于 快速 增长 和 慢 速 增长 阶段 ，a 分 别 取 0.75 和 0.42。 

瞳 曼 的 增长 分 为 早期 的 快速 增长 和 晚期 的 慢 速 增长 两 个 阶段 。 由 于 早期 的 剧烈 弛 豫 作 

用 ， 物 质 可 以 得 到 充分 的 混合 ， 位 力 半 径 (ras) 和 特征 半径 (rs) 便 可 以 同步 增长 ， 因 此 两 者 

的 比值 聚集 度 几 乎 保持 不 变 。 而 到 了 晚期 以 后 ， 物 质 难 以 再 钞 入 暗 晕 核 心 ， 只 能 逐 层 登 加 到 
外 围 ， 所 以 特征 半径 不 变 ， 位 力 半 径 却 在 增加 ， 因 而 聚集 度 剧 烈 增长 。 
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— I, 4 


=Å, 
- = = r, 4,=200/Q(z) 


0.1 


a a 


PRI 5 个 高 分 辩 率 的 暗 晕 。 各 个 子 图 中 ， 实 线 表 示 rvi:， 短 虚线 表示 raooc， 长 虚线 表示 raoom; 下 面 
蓝 色 虚线 表示 模拟 中 的 rs*， 红 色 实 线 表 示 模 型 预测 的 ro; 红色 实心 点 标记 的 转折 点 。 
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m 
A 
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RS 
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x 
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ik: Edu 

YE: 


彩色 的 线 是 随机 选取 20 Aii, AS E ae [0 o EA A E AR Ph 
Zhao ™ 提出 的 模型 所 预测 的 聚集 


图 18 ， 暗 晕 半径 和 内 部 特征 半径 的 演化 呈 
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黑色 实 线 表 示 ， 虚 线 是 根据 


[三 
Zo 


图 19 ” 主 分 支 上 的 瞳 晕 质 量 增长 历史 
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Poole 等 人 ”给 出 的 FoF 瞳 晕 和 子 结构 的 质量 随 红 移 的 演化 过 程 ， 如 图 四 所 示 。 他 们 
分 别 选 取 了 数值 模拟 中 红 移 z = 0.5,2 = 1,2 = 2 时 质量 大 于 109 Mo - h-* 的 FoF 系统 ， 
FoF 暗 晕 的 质量 随时 间 增 加 ， 子 结构 的 质量 则 先 增加 到 峰值 而 后 下 降 。 


a) b) 


YE: a), b), c) 分 别 对 应 zselect = 0.5, 1.0, 2.0 时 ， 质 量 大 于 1013 Mo - hot 的 FoF 暗 晕 的 情况 。 各 个 子 图 


c) 


上 部 的 绿色 和 蓝 色 部 分 分 别 表 示 的 是 FoF 暗 晕 和 相应 的 5 个 最 大 质量 子 晕 的 质量 演化 ， 其 中 阴影 部 分 表示 68% 
的 置信 区 间 ， 线 表示 分 布 的 峰值 概率 。 子 图 下 部 表示 的 则 是 子 景 达到 暗 晕 大 小 的 峰值 时 ( 蓝 绿色 )、 吸 积 落 入 FoF 


时 ( 品 红 色 )、 最 终 并 合 时 (GERE) 的 时 间 分 布 。 


图 20  FoF 瞳 晕 及 其 子 结构 的 质量 随时 间 的 演化 


4.2.3 静态 半径 的 情况 
Cuesta 等 人 "给 出 的 静态 半径 与 位 力 半径 之 比 随 着 时 间 (尺度 


因子 ) 的 演化 关系 ， 如 图 


DO 所 示 。 可 以 看 出 ， 以 位 力 半径 来 衡量 ， 小 质量 暗 晕 的 静态 区 域 先 增加 到 最 大 值 而 后 减 小 ， 
星系 质量 的 瞳 晕 正在 靠近 这 一 最 大 值 ， 星 系 团 大 小 的 瞳 晕 的 静态 半径 仍 在 增长 中 。 


Z 


— — 10" Me « M < 1091/1 Mo 
— — 10 Me« M < 10? 5171 Me 
----10951) Mo< M < 10917! Mo 


poe op E E A | 


1 | 1 | 1 | h 
03 04 05 06 07 0.8 0.9 
a=(l+z)? 


iE: 短 虚 线 为 星系 团 质量 大 小 的 暗 晕 ， 长 虚线 为 星系 质量 大 小 的 暗 晕 ， 实 线 为 小 质量 暗 晕 。 


图 21 ”静态 半径 与 位 力 半 径 的 比值 随时 间 的 演化 图 
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4.2.4. 回 溅 半径 的 情况 


40 卷 


O'Neil 5& AP" £i ff (1013 ~ 1015) Mc, BSE RE EL ( 见 图 22), WAT Ius- 
trisTNG Ji H^ PAY N 体 模拟 TNG300-1-DM (图 a)) 和 流体 模拟 TNG300-1 (El b), e)) 的 


结果 进行 比较 ， 后 者 给 出 了 只 用 暗物质 (图 b)) 和 所 有 物质 (图 c)) 的 情况 ， 而 前 者 


只 有 了 上 暗 物 


质 。 总 的 来 说 ，N 体 模拟 与 流体 模拟 得 到 的 Rp 的 结果 相符 合 。 另 外 ， 他 们 还 与 More 等 
人 ” (图 中 虚线 ) 以 及 Diemer” (图 中 点 线 ) 的 解析 模型 进行 了 比较 ， 虽 然 解析 模型 只 适用 于 


者 物质 ， 但 为 了 方便 参考 ， 在 各 个 分 图 中 都 列 出 。 


-一 14.5~15.0 —13.5~14.0 
-一 14.0~14.5 -一 13.0~13.5 


—— IllustrisTNG 
More 15 
-=== Diemer 20 


0.3 04 06 08 0.3 04 06 08 
Z 
a) b) c) 


0.2 04 06 0.8 


VE: a) EN 体 模拟 的 暗物质 的 情况 ，b) 是 流体 模拟 中 暗物质 的 情况 ，c) 是 流体 模拟 中 的 所 有 物质 的 情况 。 上 


半 部 分 表示 回 溅 半径 和 红 移 的 关系 ， 下 半 部 分 表示 的 是 和 N 体 模拟 的 情况 相 比 的 分 数 差 分 。 不 同 颜色 表示 不 同 


的 质量 区 间 ， 实 线 表示 数值 模拟 的 结果 ， 虚 线 和 点 线 则 是 理论 模型 的 预测 ， 线 周围 的 阴影 带 则 表示 第 16 个 和 第 


84 个 百 分 位 误差 。 


图 22 ”不 同 质量 区 间 的 暗 晕 回 溅 半径 随 红 移 的 演化 呈 


5 总 结 与 展望 


本 文 总 结 了 暗物质 量 边 界 的 几 种 定义 方式 。 根 据 球 夫 缩 模型 所 提出 的 球形 高 密 算 法 ， 将 
音 尝 定义 为 密度 达到 一 定 值 以 上 的 球形 高 密度 区 ， 常 用 的 是 位 力 半径 和 位 力 质 量 。 这 种 定义 


不 能 反映 瞳 坚 的 真实 情况 ， 会 出 现 所 谓 的 “ 伪 演 化 ”现象 ， 即 根据 球形 高 密 算 法 ， 暗 党 质 量 


和 半径 在 增长 的 过 程 中 ， 相 应 的 物理 过 程 并 没有 明显 的 变化 。 除 了 SO 算法 以 外 ， 常 用 的 瞳 


晤 定义 方法 还 包括 著名 的 FoF 算法 ， 将 数值 模拟 中 距离 小 于 一 定 值 的 粒子 连 起 来 ， 


AN Eo 


得 到 一 


后 来 人 们 为 此 进行 探索 ， 提 出 了 以 下 方法 。 根 据 瞳 坚 静 态 区 间 的 最 远 位 置 ， 定 义 静 态 
半径 ， 具 有 一 定 的 优点 。 通 过 研究 瞳 晕 的 外 围 ， 发 现 密度 轮廓 的 斜率 剧烈 下 降 ， 比 传统 的 
NFW 轮廓 和 Einasto 轮廓 更 加 陡峭 ， 据 此 提出 回 溅 半径 的 概念 。 本 文 介绍 了 两 种 算法 来 计 


4 期 RZF: 


2] A E WFR FEE 


算 回 溅 半径 ， 


到 密度 轮廓 的 斜率 最 陡峭 的 点 来 拟 合成 | 
特征 耗竭 半径 的 概念 ， 可 以 更 好 地 体现 暗 晕 的 动力 学 特 复 
径 很 好 地 反映 了 暗 晕 在 边界 处 物质 流动 情况 ， 其 意 


TS 


小 。 其 中 的 内 部 耗竭 
清晰 。 
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种 是 追踪 粒子 的 轨道 远 心 点 的 SPARTA BYE; 另 一 种 是 追踪 各 个 方向 ， 找 
| 面 的 SHELLFISH 算法 。 
E, P 


音 尝 的 内 部 耗竭 半径 和 
E [n] SR f FEA BEAK 
义 比 回 溅 半径 更 


本 文 也 介绍 了 暗 尝 边界 和 观测 相关 的 部 分 。 密 度 相 关 的 边界 比较 容易 观测 ， 在 用 SZ 效 


应 计算 星系 团 质 量 和 大 小 的 时 候 ， 常 党 


使 用 Msoo 和 相应 的 Rso00。 而 回 溅 


# 径 和 耗竭 半径 则 


相对 较 难 观测 。 但 是 当 人 们 测量 星系 团 周 围 的 星系 的 数 密度 轮廓 ， 及 用 弱 引力 透镜 测量 物质 
质量 轮廓 时 ， 也 可 以 得 到 如 同 回 溅 特征 那样 斜率 剧烈 下 降 的 特征 。 银 河 系 的 内 部 耗竭 半径 的 


位 置 ， 也 有 学 者 将 其 测 出 。 
ie LIF 


合 过 程 由 并 合 树 所 描述 ， 由 位 力 质量 和 半径 可 以 看 出 ， 其 增长 分 为 快速 增长 的 
早期 和 慢 速 增长 的 晚期 。 而 伪 演 化 便 处 在 晚期 慢 速 增长 阶段 ， 如 图 TH 所 示 。 


以 前 人 们 对 于 暗 坚 的 研究 主要 集中 在 内 部 区 域 ， 因 为 与 星系 的 形成 有 关 。 而 近 些 年 外 围 


的 研究 更 加 受到 关注 ， 这 一 点 应 当 与 回 
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The Recent Progress of the Boundaries of 
Dark Matter Haloes 


SONG Jia-yu'? 


(1. Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chi- 
nese Academy of Sciences, Shanghai 200030, China; 2. University of Chinese Academy of Sciences, 
Beijing 100049, China) 


Abstract: Dark matter haloes, whose boundary definition is controversial for a long time, 
are the building blocks of our universe. We review the explorations of this issue. Tradition- 
ally, the spherical overdensity definition, which defines a dense enough spherical region as 
a halo, can lead to the problem of pseudo-evolution. That is, according to its definition, a 
halo’s mass and radius grow many times but without much happened in its inner structure 
and physical processes. Another common method is the famous friends-of-friends algorithm, 
which links particles nearby. Both of these two algorithms are early and classical methods, 
which are fast to run but have little physical meanings. Besides, new methods can always 
been put forward. The static radius, which is defined at the outermost location of the static 
region, has been introduced. The splashback radius, is proposed because of the discovery 
that the logarithmic slope of the density profile of haloes decreases steeply at outer region. 
We then introduce two algorithms to calculate the splashback radius and show some obser- 
vational evidence of it. The inner depletion radius and characteristic depletion radius are 
defined based on halo’s dynamical properties, the former has been measured through our 
Milky Way. There’s also some connection between splashback radius and inner depletion ra- 
dius, the latter could be regarded as the outermost splashback boundary. We also introduce 
the merger tree and accretion history of dark matter haloes. Finally, we summarize different 


kinds of boundaries and show some prospects. 


Key words: dark matter halo; boundary; overdensity; splashback radius; depletion radius 


